Magnetocaloric and magnetic properties of La2NiMnO6 double perovskite
Masrour R†, , Jabar A
Laboratory of Materials, Processes, Environment and Quality, Cady Ayyed University, National School of Applied Sciences, P. O. Box 63 46000, Safi, Morocco

 

† Corresponding author. E-mail: rachidmasrour@hotmail.com

Abstract
Abstract

The magnetic effect and the magnetocaloric effect in La2NiMnO6(LNMO) double perovskite are studied using the Monte Carlo simulations. The magnetizations, specific heat values, and magnetic entropies are obtained for different exchange interactions and external magnetic fields. The adiabatic temperature is obtained. The transition temperature is deduced. The relative cooling power is established with a fixed value of exchange interaction. According to the master curve behaviors for the temperature dependence of predicted for different maximum fields, in this work it is confirmed that the paramagnetic–ferromagnetic phase transition observed for our sample is of a second order. The near room-temperature interaction and the superexchange interaction between Ni and Mn are shown to be due to the ferromagnetism of LNMO.

1. Introduction

Double perovskite of the general formula A′A″B′B″X6; where the (A′A″)-sites are occupied by rare or alkaline earth ions, (B′B″)-sites by transition metal ions and (X)-sites by oxygen or halide ions, have been known for several decades.[1,2] The structure, magnetic, and dielectric properties of the double perovskite oxides, R2NiMnO6 (R = Pr, Nd, Sm, Gd, Tb, Dy, Ho, and Y) are investigated in Ref. [3]. La2NiMnO6, exhibiting a ferromagnetic Curie temperature TC of 280 K, has received considerable attention recently,[47] not only because of the possibility of combining multiple electronic properties (ferromagnetism, magnetoresistance, magneto capacitance, and semiconductivity) in this material, but also because of the expectation that a fundamental understanding of the Ni2+–O–Mn4+ electronic interaction would provide new guidelines for designing the multiple property materials. The structural, electronic, and magnetic properties of the double perovskite Pb2FeReO6 have been studied by using the first principles.[8] The magnetic ordering was explained as the superexchange interaction between Ni2+(d8, S = 1) and Mn4+ (d3, S = 3/2) ions.[7,9] Compared with other ferromagnetic semiconductors and insulators which only exhibit magnetic ordering at very low temperatures, La2NiMnO6 demonstrates ferromagnetic ordering near room temperature, as high as 280 K.[4,10] The magnetic hysteresis displacement, thermal inversion of the magnetization, hysteresis loops jumps, and crossing branches of hysteresis loops at low magnetic fields are reviewed.[11] The electronic and magnetic structures for the ferromagnetic and antiferromagnetic states of La2MnNiO6 and the site-projected density of states and electronic dispersion curves are analyzed in Ref. [12]. La2NiMnO6 is a ferromagnetic insulator with a positive superexchange interaction between Ni and Mn cations.[5,6,1315] The influences of bismuth on the magnetic and electrical properties of La2MnNiO6 are studied in Ref. [16]. The Monte Carlo data by finite-size scaling, and the location of a line of critical points of the Ising universality class are analyzed.[17] The magneto-caloric effect of a Gd55Co25Al18Sn2 bulk metallic glass is investigated.[18] The crystal field and external magnetic field-split ground state of Dy3+ in Dy3Al5O12 are calculated based on the quantum theory.[19] In this paper, theoretical work on the dependence of the magnetization on temperature for La2MnNiO6 compound is simulated to predict magnetocaloric properties such as magnetic entropy change, relative cooling power, and heat capacity change. The thermal magnetization, dM/dT, specific heat, and magnetic entropy are given for a fixed magnetic parameter. The transition temperature TC values of La2NiMnO6 double perovskite are deduced for different values of external magnetic field. The thermal total magnetization, specific heat, and magnetic entropy are established with the external magnetic field h. We also give the adiabatic temperatures with different values of external magnetic field. Finally, the field dependence of relative cooling power (RCP) is obtained.

2. Model and formulation

The Hamiltonian of an Ising model with the La2NiMnO6 double perovskite includes nearest neighbor interactions and external magnetic field and is given as

where ⟨i,j⟩ denote the first and second nearest neighbor sites i and j; h is the external magnetic field; Jij(Ni−Ni), Jij(Mn−Mn), Jij(Mn−Ni) are the first (J1(Ni−Ni), J1(Mn−Mn), J1(Mn−Ni)) and second J2(Ni−Ni), J2(Mn−Mn), J2(Mn−Ni) nearest neighbor between the NiNi, MnMn, and MnNi respectively. The spin moments of Mn4+ and Ni2+ ions are S = 3/2 and σ = 1, respectively. We use the mean field theory to calculate the exchange interactions between Ni–Ni, Mn–Mn, and Mn–Ni. The obtained values are J1(Ni−Ni) = 22.0 K, J1(Mn−Mn) = 24 K, J2(Ni−Ni) = 19.0 K, J2(Mn−Mn) = 23.0 K, J1(Mn−Ni) = 15.0 K, and J2(Mn−Ni) = 10.0 K. The values of spin moment S and σ take ±3/2; ±1/2, and ±1,0, respectively.

3. Monte Carlo simulations

The La2NiMnO6 double perovskite as given in Fig. 1 is assumed to reside in the unit cell and the system consists of the total number of spins N = NMn + NNi, with NMn = NNi = 1024.

Fig. 1. Double perovskite La2NiMnO6 structure.

We use a standard sampling method to simulate the Hamiltonian given by Eq. (1). Cyclic boundary conditions on the lattice are imposed and the configurations are generated by sequentially traversing the lattice and making single-spin flip attempts. The flips are accepted or rejected according to a heat-bath algorithm under the Metropolis approximation. Our data are generated with 105 Monte Carlo steps per spin, discarding the first 104 Monte Carlo simulations. Starting from different initial conditions, we perform the average of each parameter and estimate the Monte Carlo simulations, averaging over many initial conditions. Our program calculates the following parameters:

4. Results and discussion

In a previous work, we have studied the magnetocaloric effect in NdMnO3 perovskite by using the Monte Carlo simulations.[20] La2NiMnO6 is a double perovskite structure and is so named because the unit cell is twice that of perovskite. They do not have the same structures. In NdMnO3, we consider the first nearest neighbor but in the La2NiMnO6 system we consider the first and the second exchange interactions for each atom.

Shown in Fig. 2(a) are the thermal total magnetizations of La2NiMnO6 and thermal magnetization of Ni and Mn. The dM/dT, specific heat, and magnetic entropy of La2NiMnO6 are given in Figs. 2(b)2(d), respectively under h = 0. The ferromagnetic (FM) transition temperature is calculated from the minimum position of the (dM/dT) versus temperature curve. The observed ferromagnetism could be explained within the framework of FM interaction between the Mn and Ni ions.[6] The transition temperature obtained is TC = 275 K. The transition from a paramagnetic (PM) state to an FM state around TC ≈ 275 K for the parent phase as reported in Refs. [6] and [16]. The near room temperature interaction and superexchange interaction between Ni+2 and Mn+4 are shown to be due to the ferromagnetism of LNMO. This behavior is observed in the previous work.[21] The double perovskite La2NiMnO6 is ferromagnetic at transition temperature TC = 280 K, which has been predicted according to the Goodenough–Kanamori rules and confirmed by the experimental data.[15,2224]

Fig. 2. Thermal magnetizations (a), dM/dT (b), the specific heat (c), and magnetic entropy (d) of La2NiMnO6 for h = 0.

The magnetic entropy increases with increasing temperature until the transition temperature TC, which corresponds to the maximum value of the magnetic entropy. The thermal total magnetization and dM/dT of La2NiMnO6 are shown in Figs. 3(a) and 3(b), respectively for external magnetic field h = 1, 4, 7, and 10 T. The total magnetization increases with the external magnetic field increasing. On the other hand, Blasse argued that ferromagnetism in La2NiMnO6 is entirely due to Mn4+–O–Ni2+ superexchange interactions,[16] which was later concluded[25] from the measurement of 55Mn nuclear magnetic resonance line width. The curves of magnetic entropy versus the temperature of La2NiMnO6 are given in Fig. 4 for different values of external magnetic field h = 1, 4, 7, and 10 T. The magnetic entropy increases with increasing temperature until the transition temperature TC and becomes constant for different values of external magnetic field.

Fig. 3. Thermal total magnetizations (a) and dM/dT curves (b) of La2NiMnO6 for external magnetic field h = 1, 4, 7, and 10 T.
Fig. 4. Curves of magnetic entropy versus the temperature of La2NiMnO6 for different values of external magnetic field h = 1, 4, 7, and 10 T.

We show in Fig. 5(a) the temperature dependences of the magnetic entropy for different values of external magnetic field h = 1, 4, 7, and 10 T. The value of ΔSm is negative in the entire temperature range and is extended over a wide range of temperature around the Curie temperature, which is useful for above room temperature magnetic refrigeration. The values of −ΔSm reach 8.44, 15.45, 20.65, and 25 J/(kg·K) at TC ≈ 275 K when the applied magnetic fields are 1, 4, 7, and 10 T, respectively. The paramagnetic–ferromagnetic transition is extremely sharp, indicating its second-order phase character as given in Fig. 5(b). One can see that all the ΔSm(T,h) data points are collapsed into a universal curve in the whole temperature range as given in Ref. [26]. As is well known, the magnetic order is a second order phase transition in the La2NiMnO6. Therefore, the corresponding lattice and electronic entropy change will be much smaller. These also rule out the possibility that the magnetic entropy change is partially compensated for by the lattice and electronic contributions. This behavior is observed in previous work.[15]

Fig. 5. Temperature dependences of the magnetic entropy change (a) and master curve behaviors of the curve each as a function of temperature (b) for La2NiMnO6 for different values of external magnetic field h = 1, 4, 7, and 10 T.

A large change in magnetic entropy due to adiabatic magnetization is observed, which argues well for potential applications as magnetic refrigerant around liquid nitrogen temperature. Luo et al. suggested that the large magnetic entropy change in perovskite manganites could originate from the spin–lattice coupling in the magnetic ordering process.[27] We have also presented in Fig. 6, the temperature dependences of the adiabatic temperature in La2NiMnO6 for different values of external magnetic field h = 1, 4, 7, and 10 T. A large magnetic entropy change is observed and the changes of magnetic entropy and adiabatic temperature at moderate magnetic fields are quite significant. A maximum adiabatic temperature change about 275 K is obtained to be close to the transition temperature with applied field h = 1, 4, 7, and 10 T. These results indicate that our samples are very good substances for magnetic refrigeration applications. To explain a large magnetic entropy change in perovskite manganite, the Zener double exchange model has been strongly recommended.[28,29] A maximum entropy, −ΔSM = 20.54 J · kg−1 · K−1 and adiabatic temperature change 4.84 K are comparable to these values: 35.5 J · kg−1 · K−1 (∼ 24 J · kg−1 · K−1, and 10.5 K (6.5 K) is observed in Gd2NiMnO6 (Gd2CoMnO6) for a magnetic field change of 7 T at low temperatures.[30]

Fig. 6. Temperature dependences of the adiabatic temperature change for La2NiMnO6 for different values of external magnetic field h = 1, 4, 7, and 10 T.

Figure 7 shows the field dependence of relative cooling power (RCP) on external magnetic field for La2NiMnO6. It can be seen from this figure that RCP increases monotonically as the field increases, and reaches the value of 905 J·kg for h = 10 T. The RCP is another important factor for assessing the usefulness of a magnetic refrigerant material.[31] The RCP depends not only on the magnitude of ΔSm, but also on the temperature dependence of ΔSm (e.g., the full width at half maximum of the ).[32] The values obtained around TC = 275 K increase with applied field increasing. All the parameters related to the magnetocaloric effect, such as the magnetic entropy change, relative cooling power, and refrigerant capacity depend strongly on the magnitude of the applied magnetic field h and exchange interaction between Mn and Ni ions.

Fig. 7. Field dependence of relative cooling power (RCP) for La2NiMnO6.

Finally, we show in Fig. 8 the hysteresis loops of La2NiMnO6 for different temperatures T = 260, 275, and 285 K. The coercive field decreases with temperature increasing. The same behavior has been observed experimentally.[33] The system exhibits the superparamagnetic behavior for TTC. However, although the large abrupt change in magnetization causes correspondingly a giant magnetic entropy change, this appears at the expense of thermal and magnetic hysteresis, which should be avoided in order to be able to apply these materials in refrigerators appliances. The coercive field in this double perovskite is sample-dependent, which might indicate that the role of structural defects such as antiphase boundaries[34] acting as domain pinning centres can be very important in such highly anisotropic compounds. Antiphase boundaries in double perovskite can occur when there appears a lattice shift equal to half a lattice parameter along a crystallite plane. This gives rise to a plane formed by antiferromagnetic Ni–O–Ni bonds or Mn–O–Mn bonds. The extremely large magnetic entropy change may be attributed to a large number of weakly interacting spins. These indicate that the material is a promising magnetocaloric effect (MCE) candidate for low temperature application, and possibly could make ultra-low temperatures easily achievable for most laboratories and for space application as well.

Fig. 8. Hysteresis loops of La2NiMnO6 for different temperatures T = 260, 275, and 285 K.
5. Conclusions

The magnetocaloric effect and the critical behavior of the double perovskitesLa2NiMnO6 are investigated by Monte Carlo simulations. The maximum of the magnetocaloric effect corresponds to transition temperature TC. The La2NiMnO6 shows the second order PM–FM phase transition around 12 with a large magnetic entropy change over a wide range of temperature, which is comparable to the magnetic entropy changes of other manganites. The changes in magnetic entropy and the adiabatic temperature are also significant at moderate magnetic fields. Therefore, a great deal of attention has been focused on finding new materials with a large MCE and small magnetic hysteresis.[35,36] These results are shown in our material. The exchange interaction between Mn–Ni atoms is the main mechanism that controls the magnetic and transport properties of manganite and are chosen the positive values as given in Refs. [5], [6], [13], and [15]. Finally, in order to study the nature of the paramagnetic–ferromagnetic phase transition, a master curve behavior for the temperature dependence of predicted for different maximum fields is proposed and the hysteresis loop is obtained for the temperatures superior and inferior to those of the transition temperature.

Reference
1Flerov I NGorev M VAleksandrov K STressaud AGrannec JCouzi M 1998 Mater. Sci. Eng. 24 81
2Anderson M TGreenwood K BTaylor G APoeppelmeier K R 1993 Prog. Solid State Chem. 22 197
3Booth R JFillman RWhitaker HNag ATiwari R MRamanujachary K VGopalakrishnan JLofland S E 2009 Mater. Res. Bull. 44 1559
4Rogado N SLi JSleight A WSubramanian M A 2005 Adv. Mater. 17 2225
5Dass R IYan J QGoodenough J B 2003 Phys. Rev. 68 064415
6Joseph V LJoy P ADate S KGopinath C S 2002 Phys. Rev. 65 184416
7Sanchez M CGarcia JBlasco JSubias GPerez-Cacho J 2002 Phys. Rev. 65 144409
8Zhang YJi VXu K W 2012 Physica 407 2617
9Wang TShi W ZFang X DDong W WTao R H 2010 J. Sol–Gel Sci. Technol. 53 655
10Hashisaka MKan BMasuno ATakano MShimakawa Y 2006 Appl. Phys. Lett. 89 032504
11Moure CPena O 2013 J. Magn. Magn. Mater. 337�?38 1
12Matar S FSubramanian M AVillesuzanne AEyert VWhangbo M H 2007 J. Magn. Magn. Mater. 308 116
13Goodenough J BWold AArnott R JMenyuk N 1961 Phys. Rev. 124 373
14Blasse G 1965 J. Phys. Chem. Solids 26 1969
15Choudhury DMandal PMathieu RHazarika ARajan SSundaresan AWaghmare U VKnut RKaris ONordblad PSarma D D 2012 Phys. Rev. Lett. 108 127201
16Nautiyal PMotin Seikh MPralong VKundu A K 2013 J. Magn. Magn. Mater. 347 111
17Lv J PDeng Y JChen Q H 2011 Phys. Rev. 84 021125
18Ding DZhang Y QXia L 2015 Chin. Phys. Lett. 32 106101
19Xue DYing Z GSuo X WJia H YFeng H 2015 Acta Phys. Sin. 64 177502 (in Chinese)
20Masrour RJabar ABenyoussef AHamedoun MHlil E K 2016 J. Magn. Magn. Mater. 401 91
21Zhou SShi LYang HZhao J 2007 Appl. Phys. Lett. 91 172505
22Wang XSui YLi YLi X LWang Z YLiu ZSu WTang J 2009 Appl. Phys. Lett. 95 252502
23Wold AArnott R JGoodenough J 1958 Appl. Phys. Berl. 29 387
24Goodenough J B 1955 Phys. Rev. 100 564
25Sonobe MAsai K 1992 J. Phys. Soc. Jpn. 61 4193
26Manh T VHo T AThanh T DPhan T LPhan M HYu S C 2015 IEEE Trans. Magnetic. 51 2400304
27Luo XSun Y PWang BZhu X BSong W HYang Z RDai J M 2009 Solid State Commun. 149 810
28Guo Z BDu Y WZhu J SHuang HDing W PFeng D 1997 Phys. Rev. Lett. 78 1142
29Wang Z MNi GXu Q YSang HDu Y W 2001 J. Appl. Phys. 90 5689
30Murthy K JChandrasekhar K DMahana STopwal DVenimadhav A 2015 J. Phys. D: Appl. Phys. 48 355001
31Tang TGu K MCao Q QWang D HZang S YDu Y W 2000 J. Magn. Magn. Mater. 222 110
32Kitanovski AEgolf P W 2006 Int. J. Refrig. 29 3
33Guo H ZBurgess JAda EStreet SGupta AIliev M NKellock A JMagen CVarela MPennycook S J 2008 Phys. Rev. 77 174423
34Yu XAsaka TTomioka YTsuruta CNaai TKimoto KKaneko YTokura YMatsui Y 2005 J. Electron Microsc. 54 61
35Cubrovic MZaanen JSchalm K 2009 Science 325 439
36Pfleiderer CMcMullan G JJulian S RLonzarich G G 1997 Phys. Rev. 55 8330